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Abstract 

J^ I We deploy numerical semidefinite programming and conversion to exact rational in- 

ly^ ■ equalities to certify that for a positive semidefinite input polynomial or rational func- 

^N , tion, any representation as a fraction of sums-of-squares of polynomials with real co- 

. I efficients must contain polynomials in the denominator of degree no less than a given 

^12 ' input lower bound. By Artin's solution to Hilbert's 17th problems, such representations 

p^ ■ always exist for some denominator degree. Our certificates of infeasibility are based on 

the generalization of Farkas's Lemma to semidefinite programming. 

The literature has many famous examples of impossibility of SOS representability 
k>( ■ including Motzkin's, Robinson's, Choi's and Lam's polynomials, and Reznick's lower 

;_( , degree bounds on uniform denominators, e.g., powers of the sum-of-squares of each 



variable. Our work on exact certificates for positive semidefiniteness allows for non- 
uniform denominators, which can have lower degree and are often easier to convert 
to exact identities. Here we demonstrate our algorithm by computing certificates of 
impossibilities for an arbitrary sum-of-squares denominator of degree 2 and 4 for some 
symmetric sextics in 4 and 5 variables, respectively. We can also certify impossibility 
of base polynomials in the denominator of restricted term structure, for instance as in 
Landau's reduction by one less variable. 
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1. Introduction 

The Farkas Lemma of linear programming can be employed to construct certificates of infea- 
sibility, in its simplest form of an inconsistent system of linear equations [Giesbrecht, Lobo, 
and Saunders 1998], in linear programming of a system of linear inequalities, and in semidef- 
inite programming of a system of linear equations with semidefiniteness constraints on the 
solution. A polynomial is not a sum-of-squares of polynomials (SOS) if the corresponding 
semidefinite program is infeasible. Thus the Farkas Lemma produces a certificate that a 
polynomial is not an SOS, the separating hyperplane [Ahmadi and Parrilo 2011]. 

Motivated by our SOS certificates for global optima of polynomials and rational functions 
[Kaltofen, Li, Yang, and Zhi 2012; Kaltofen, Yang, and Zhi 2009; Hutton, Kaltofen, and 
Zhi 2010] (see also [Powers and Wormann 1998; Harrison 2007; Peyrl and Parrilo 2008] 
for earlier work), we extend those impossibility certificates to Hilbert-Artin representations 
of a given denominator degree: by Emil Artin's Theorem [Artin 1927], every real positive 
semidefinite rational function is a fraction of two sums-of-squares of polynomials. We write 
for an f{Xi, . . . , X„) G K{Xi, . . . , X„), where i^ ^ Q is a subfield of the real numbers, 

fho a vei,...,eneM:/(ei,...,en)^o. 

Note that at a real root of the denominator of / its value is undefined, hence ^ 0. Artin's 
original theorem stipulates that 

1 ' 
V/^0: 3m,...,Ui,weK[X,,...,X„]: f = —y2ul (1) 

i=l 

If / is a polynomial G K[Xi, . . . , X„], one may eliminate one variable from the denominator 
w, that is, construct from (1) a representation with tfnew ^ K[Xi, . . . ,X„,-i], which Artin 
in his 1927 paper attributes to Edmund Landau. The reduction can be accomplished with 
the same number of /new = / squares, which in [Rajwade 1993] is attributed to J. W. S. 
Cassels. In both constructions the degree of Wnew is substantially larger than that of w. As is 
customary, if a positive semidefinite polynomial / allows a representation (1) with w = 1, we 
shall call / a sum-of-squares (SOS). In general, however, as already David Hilbert has shown 
in 1888, positive semidefinite polynomials are not SOS [Chesi 2007; Blekherman 2009]. 
In order to minimize the numerator and denominator degrees, we seek 



2^i=l 



ul 



Ui,...,uu Vi,. . .,vi' G M[Xi, . . . , X„] such that / = —j, "-. (2) 



'' „2' 



E,=i^, 



J 



We shall call (2) a Hilbert-Artin representation of /, which constitutes an SOS proof for 
/ ^ 0. By allowing an SOS as the denominator polynomial, one then can construct such 
proofs with a possibly smaller degree than the common denominator w"^ in (1). For instance, 
for the Motzkin polynomial maxj{deg(wj)} < 1 suffices in (2), but deg{w) < 1 is impossible 
in (1) [Kaltofen, Li, Yang, and Zhi 2012, Section 1]. 

It is not known if minimal degree denominator SOSes can always have coefficients in 
K, as is the case in Artin's original theorem (1). A special case is when / is an SOS of 
polynomials {w = 1), and the existence of Ui with all coefficients in K for all i is conjectured 



(Sturmfels; cf. [Hillar 2009; Kaltofen 2009; Quarez 2009; Scheiderer 2009]). Our method can 
certify "absolute" impossibility by SOSes, that is, for coefficients from all possible subfields 
of M. Our certificates are rational, that is, they have their scalars in Q. The problem whether 
there exists a representation of a given degree with coefficients in Q appears to be decidable 
[Safey El Din and Zhi 2010]. 

As in [Kaltofen, Li, Yang, and Zhi 2012], we compute our certificates, the separating 
hyperplanes in Farkas's Lemma, by first computing a numerical approximation numerical 
semidefinite program solver and then converting the numerical scalars to exact rational num- 
bers. For ill-posed polynomials (see Example 5.4 below), high-accuracy semidefinite program 
solver [Guo 2009] is needed. The separating hyperplane is the strictly feasible solution to 
a semidefinite program whose objective function tends to — oo. We compute such a strictly 
feasible solution by the Big-M method [Vandenberghe and Boyd 1996]. The semidefinite pro- 
grams in [Ahmadi and Parrilo 2011] and ours certify infeasibility of SOSes, which has been 
generalized to infeasibility of arbitrary linear matrix inequalities [Klep and Schweighofer 
2011]. 

We have tested our method on polynomials from the literature. In particular, we show that 
the SOS proofs of positive semidefiniteness in [Kaltofen, Yang, and Zhi 2009] indeed require 
denominators for three polynomials. The ArtinProver program [Kaltofen, Li, Yang, and Zhi 
2012] successfully introduced denominators not only for purpose of handling inequalities that 
do not allow a polynomial SOS proof, but also for avoiding possible non-rational SOSes, 
to which the semidefinite program solvers may have converged in the case where the Gram 
matrix is intrinsically rank deficient (our "hard case" [Kaltofen, Li, Yang, and Zhi 2012]). 
Our impossibility certificates show that for the proof of the Monotone Column Permanent 
conjecture in dimension 4, actually the former is the case. 

A final problem is to explicitly construct a positive semidefinite polynomial for which 
the Hilbert-Artin representation (2) must have deg(X]j v'j) > 4. Bruce Reznick in 2009 has 
kindly provided us with the challenges raised in [Choi, Lam, and Reznick 1987, Section 7]: 
how necessarily high must be the powers {xf + ■ ■ ■ +x'^Y in the (uniform) denominators such 
that a family fn,k (see Example 5.3 below) of even symmetric sextics in n variables is an 
SOS, where 2 < k < n — 27 In [Choi, Lam, and Reznick 1987] it is proven that for /4 2 one 
has r = 2. We can compute certificates that show that for /42, /5,25 f6,2, the degree lower 
bound > 4 and for /ss, /e4, the lower bound > 6 even hold for any denominator ^ -v"^ in 
(2). 

Notation Throughout this paper, N denotes the set of nonnegative integers and we set 
NJ^ = {a G N" I |a| = Er=i "* < for ^ ^ N. R[X] = M[Xi,...,X„] denotes the ring 
of polynomials in variables X = [Xi, . . . ,Xn) with real coefficients. Given a polynomial 
/ = Ea/«^r' • • -^n" e R[X], let supp(/) = {Xfi . ..X^" I ca ^ 0}, i.e., the set of the 
support terms of /. Denote by deg(/) the total degree of /. Given n > 1 and e > 0, let 
Terms[X; deg<e] = {^i^ . . . X^" \ Yl'i=i ^i ^ ^j, i.e., the set of all terms of total degree < e 
in the n variables X. For a given subset T C Terms[X; deg<e], we introduce the following 
notation for a term-restricted SOS, SOSr = {^"^1 I "^j ^ M[X], supp(fj) C T}, and the 
following notation for a denominator term-restricted Hilbert-Artin representation, 

SOS/SOSr ={Y1 "i / Yl ^i I "" ^J ^ ^[X],yj ■■ supp(?;j) C t\. (3) 



Finally, we write the shorthand S0S/S0Sdeg<2e = SOS/SOSTerms[X;deg<e]- 

By SM*^^*^ we denote the subspace of real symmetric k x k matrices. For a matrix W e 
gj^fcxfc^ lY y Q means W is positive semi definite. The bold number zero denotes the zero 
matrix, and / denotes the identity matrix. 

2. Hilbert-Artin representation of positive semidefinite 
polynomials 

2.1. Rational Function Sum-of- Squares and Semidefinite Program- 
ming 

For a given subset T C Terms[X; deg<e], note that / G SOS/SOSr if and only if 

I I' 

* j 

for some polynomials Ui{X),Vj{X) G M.[X] with supp(fj) G T. Consider the following set: 

mTcvms[X;dcg<d]'^W^^^mTerms[X;dcg<d] = f {X) ■ TUr'^W^'^^mr 



[W^'\W^ 



lyw y 0, iy[2i y o, Tt{w^^^) = i r > (4) 



where m-j- and ?TiTcrms[X;deg<d] denote the column vectors which consist of the elements in 
T and Terms[X; deg<(i], respectively. Here and hereafter, we let d = \e + deg(/)/2], and 
therefore, 

|^a+/3 I x^'^X^ G Terms[X;deg<d]} D {X"+^+^ | X^ G supp(/),X",X^ G T} . 
The last constraint Tr(14^[^l) = 1 is added to enforce that W^^^ 7^ 0. 

Proposition 2.1 We have f ^ SOS/SOSr if ^.^d only if the set (4) is empty. 

Now we review the following standard Semidefinite Program, (SDP) (see [Vandenberghe 
and Boyd 1996]), 



(5) 



sup —C • W inf b^y 

s.t. A*W = h, t = l---l, s.t. C + E!=i2/iA^0. 

Who. 

For symmetric matrices C, W, the scalar product in M"^" space is defined as 

n n 

C*W ={C,W) = ^^QjWij = TiCW. 
i=l j=l 

Let 

mTerms[X;dcg<d]^W^'^mTcrms[X;dcg<d] = ^(G^"! • W^'^)X'', (6) 



where Gf°l are scalar symmetric matrices and X°' are all possible terms appearing in the 
polynomial of degree < 2d. Similarly, let 

where H^^"^ are symmetric matrices and X^ are all possible terms appearing in the product 
of (— /(X)) with a polynomial of degree < 2e. Now we consider the following block SDP: 



sup 

VygSRfcxfc 



S.t. 



c*w 



A\-' 



W 



A*W 



H^^O. 



(7) 



where k = r/) + (T) 



C :-- 











w 






^w 



Gt' 



i/I' 



A 







and a ranges over Ng^. The matrix C can be chosen as a random symmetric matrix. We set 
it to be a zero matrix only for the convenience of discussions below. For all block positive 
semidefinite matrices appearing in the present paper, we use the symbol * to indicate that 
the associated elements could be any real numbers such that the whole matrices are still 
positive semidefinite and leave some positions blank to indicate that the associated blocks 
are zero matrices. 

Proposition 2.2 We have f ^ SOS/SOSr if and only if SDP (7) is infeasible. 



2.2. Dual Problem and Certification 

Before we consider the dual problem of (7), let us review some definitions about moment 
matrices and localizing moment matrices. Given a sequence y = {ya)aG'M" ^ M^", its moment 
matrix is the (infinite) real symmetric matrix M{y) indexed by N", with (a, /3)th entry ya+i3, 
for a,/3 G N". Given an integer t > 1 and a truncated sequence y = {ya)a£Nri^ ^ M.^^*, its 
moment matrix of order t is the matrix Mt{y) with {a, /3)th entry ya+i3, for a, /3 G N". For a 
given polynomial q G M[X], if the (z, j)th entry of Mt{y) is y^, then the tth localizing moment 
matrix of q is defined by 

Mt{qy){i,j) := ^qaya+i3- 

a 

More details about moment matrices, see [Lasserre 2001, 2009; Laurent 2009]. 
According to (5), the dual problem of (7) is 



inf s 

(j/,s)GlR™+l 



S.t. M{y,s) hO, 



(8) 



where 

^^'^- [ M,{{-f)y) + sI 

y '■= {ya)a<^N" £ ]R^2d and m = ("^^ ), Md{y) is a truncated moment matrix of order d and 
Me{{— f)y) is the eth locahzing moment matrix. 

The next lemma shows that the problem (8) is strictly feasible. The proof is similar to 
the one given in [Lasserre 2001, Propositions. 1]. 

Lemma 2.3 There exists {y, s) G ]R'"+^ such that Mdiy) >- and M^i^-fy) + 5/^0. 

Proof. Let /x be a probability measure on M"' with a strictly positive density h with respect 
to Lebesgue measure such that 

ya ■= / X"dfi < oo. 



For any polynomial q{X) G ]R[X] with supp(g) G Terms[X; deg<(i], let vec(g) denote its 
sequence of coefficients in the monomial basis Terms[X; deg<d]. We have 



(vec(g)^',Mrf(y)vec(g)) = / q{xf^i{dx) 

q{xYh{x)dx 



> whenever q ^ 0, 

which implies Md{y) >- 0. Take 

5 > -Xmin{M,{{-f)y)), 

then Me{-fy) + si y 0. D 

For standard SDPs in (5), we have the following important duality fact. 

Lemma 2.4 [Alizadeh 1995, Lemma 2.3; Semidefinite Farkas Lemma] 

Let Ai G SM'""^'^ for alii = 1, . . . , / and let 6 G M'. Suppose there exists a vector y eM! such 

that X]i=i yi^i '^ 0- Then exactly one of the following is true: 

1. There exists a positive semidefinite symmetric matrix W G SM'^^'^, W ^ 0, such that 
Ai • W = bi for all i = 1, . . . , I; 

2. There exists a vector y eW such that X]j=i IJi^i ^ o.'i^d b^y < 0. 

We call the vector y Farkas 's certificate vector of infeasibility. For other forms of the Farkas 
Lemma, see [Dattorro 2011, Section 4.2]. 

Note that Lemma 2.3 implies that the assumption in the Farkas Lemma 2.4 is satisfied 
in SDPs (7) and (8). Then we have our main result: 

Theorem 2.5 Given a polynomial f G Q[X] and an integer e > 0, let d = [e + deg(/)/2], 
then for any subset T C Terms[X; deg<e], the following are equivalent: 

I. f i SOS/SOSr, 



2. There exists a rational vector y = {i/a) G Q"* with m = ("^J ) such that Mdijj) h and 
Meify) -< 0. 

Proof. By employing the Farkas Lemma 2.4 to SDPs (7) and (8), we have that / ^ 
SOS/SOSr if and only if there exists p' = {y', s') E M™+^ for (8) such that M{y', s') >z and 
s' < 0. Now we prove that p' can be chosen to be rational. 

Let p = {y,s) be the strictly feasible point constructed in Lemma 2.3. For < t < 1, 
let 1/ = (1 — t)y' + ty and s = (1 — t)s' + ts, then M{y, s) >- 0. Since s' < 0, it is always 
possible to choose a rational number t such that s < 0. Then there exists e > such 
that for all p = {y, s) G Bp{e) where Bp{e) is a ball with center p and radius e, we have 
M{y,s) y 0. Taking e < ^\s\, there always exists a point p = {y,s) G Bp{e) such that 
p G Q"'+\ Mdijj) h 0, Mei-fy) + sJ ^ and s < which implies Meify) ^0. D 

2.3. Moment matrices and linear forms on M[X] 

In this section, we give an interpretation of our infeasibility certification using linear forms 
onM[X]. 

Given y G M^", we define the linear form Ly G (]R[X])* by 

Lyif) := y^Yecif) = J^Vafa for / = J^/.X" G R[X], (9) 

a a 

where vec(/) denotes its sequence of coefficients. 

Lemma 2.6 [Laurent 2009, Lemma 4.1] Let y G M.^" , Ly G (]R[X])* the associated linear 
form, and let f,g,hE M.[X]. 

1- Ly{fg) = vec{f)^M{y)vec{g); m particular, Ly{f^) = vec{f)^M{y)vec{f), Ly{f) = 
vec{l)^M{y)vec{f). 

2. Lyifgh) = vec{fYM{y)vec{gh) = vec{fgfM{y)vec{h) = vec{fYM{hy)vec{g). 

Now we have the following statement which is equivalent to Theorem 2.5: 

Theorem 2.7 Given a polynomial f G Q[X] and an integer e > 0, let d = [e + deg(/)/2'l, 
then for any subset T C Terms[X; deg<e], the following are equivalent: 

1. f i SOS/SOSr, 

2. There exists a rational vector y G Q™ with m = ("^^ ) , and the associated linear form 
Ly G (]R[X]2d)* such that for any polynomials v , u G M.[X] with supp{v) G T and suppiu) G 
Terms [X; deg<(i], we have Lyifv"^) < and Lyiv^) > 0. 

Proof. By (8) and Theorem 2.5, we have that / ^ SOS/SOS7- if and only if there exists 
y G Q™ such that Md{y) ^ and Me{fy) -< 0. According to Lemma 2.6, the conclusion 
follows. n 

Now one has a better understanding that the existence of a certificate y in Theorem 2.5 
implies / ^ SOS/SOS7-. In fact if / = J2'^'i/J2'^] with supp(uj) G Terms [X ; deg<(i] and 
supp(wj) G T, then < Liji^uf) = J^^yif"^']) < which is a contradiction. 

7 



Remark 2.1 One special case is e = 0, i.e. we certify that / can not be written as a rational 
SOS. According to Theorem 2.7, / is not an SOS if and only if there is ^ G Q"' and the 
associated linear form Ly, such that Vm G M.[X] with supp(m) G Terms[X; deg<[deg(/)/2] ], 
Ly{u^) > and Ly{f) < 0. This special case has also been studied in [Ahmadi and Parrilo 
2011], in which y is referred as the separating hyperplane. 



3. Computational aspects of the certification 

3.1. Finding y by Big-M method 

Given a polynomial / G Q[X] and an integer e > 0, note that / ^ SOS/SOS7- if and only if 
(7) is infeasible. From the proof of Theorem 2.5, we have 



Lemma 3.1 / ^ SOS/SOS7- if and only if s* 



-00 in 



To find a certificate y in Theorem 2.5, we need find a feasible point of the dual problem 
(8) at which the value of its objective function s is negative. We employ the Big-M method 
[Vandenberghe and Boyd 1996] to (7) and (8), and solve the following two modified SDPs 



r := sup 



s.t. 



C»{W 



w) 



Mw 



A[°] miW -w) 



_ A»{W 



w] 





1 



ly ^0,m; > 0, 



(10) 



inf s 

(3/,s)GR'"+i 



S.t. M{y,s) h 0, 
TrM(y,s) <M, 



:iii 



where matrices C, A^°'\ A, M(y, s) in (10) and (11) are defined as in (7) and (8). Note that 
any feasible point of (11) is also feasible to (8). As shown in [Vandenberghe and Boyd 1996], 
(10) and (11) are always strictly feasible and r* = s* ^ —00 as A^ — )■ 00. Hence a certificate 
y is obtained by solving (10) and (11) using interior-point methods. 

Algorithm 3.1 

Input: / G Q[X], e G Z>o and a subset T C Terms [X; deg<e]. 

Output: If f i SOS/SO~Sr, return a certificate y G Q"^. 

I. Reduce the problem to SDPs (7) and (8). 

II. Fix abig M eZ and modify (7), (8) to (10), (11). 

III. Solve (10) and (11) by interior-point methods until a solution pk = {y^,s^) with s^ < 
is obtained. 



IV. Find a strictly feasible point p = [y, s) of (8). 
V. Fix < t < 1 and p = (1 — t)pk + tp= {y, s) such that s < 0. 
VI. Choose a rational point p = {y, s) G B^{j)) where e < \\s\. 

Remark 3.1 In Step III, provided that the problem is of large size and not ill-conditioned, 
we can solve (10) and (11) using SDP solvers in Matlab like SeDuMi [Sturm 1999] which is 
very efficient. If the problem has small size and an accurate solution is needed. Maple package 
SDPTools [Guo 2009] is a better choice. SDPTools, in which the above algorithm has been 
implemented, is a high precision SDP solver based on the potential reduction method in 
[Vandenberghe and Boyd 1996]. 

Remark 3.2 In practice, if (10) and (11) in Step III are precisely computed by interior- 
point methods, then the floating-point solution [y^, s'^) is a highly accurate approximation of 
a strictly feasible point of (8). Hence, without Step IV, V, VI, one can expect that an exact 
certiflcate can be obtained by simply rounding (y^, s^) to a rational feasible solution to (8). 

3.2. Exploiting the sparsity 

To reduce computation cost, we can replace mTerms[X;dcg<d] in (4), i.e. the vector of all terms 
with degree < rf by a sparse vector containing part of mTerms[X;deg<d] due to the following 
theorem: 

Theorem 3.2 [Reznick 1978, Theorem 1] For a polynomial p{x) = XloPa^"; ^^ define C{p) 
as the convex hull of {a\ pa 7^ 0}, then we have C{p^) = 2C{p); for any positive semidefinite 
polynomials f and g, C{f) C C{f + g); if f = Y^.g] then C{gj) C |C(/). 

Definition 1 Given a polynomial f G M[X], an integer e > and a subset T C Terms[X; deg<e] 
let Cf^r be the convex hull 0/ {a G N" | a = /3 + 71 + 72, X^ G supp{f),X"'\X"'^ eT}. We 
define Qj^j- = {X°'\ 2a G C/^r}- UT = Terms[X; deg<e], we write the shorthand Gf,dog<e- 

By Theorem 3.2, we have mg^.^ C mTerms[X;deg<d] and that / G SOS/SOSr if and only if 

= mg^^^W^^^mg^.^ + (-/(X)) ■ mr^W^^^mr. 

Thus the sizes of the SDPs (7) and (8) decrease. We show a sparse version of Theorem 2.7 
below. 

Corollary 3.3 Given a polynomial f G Q[X] and an integer e > 0, let d = \e + deg(/)/2], 
then for any subset T C Terms[X; deg<e], the following are equivalent: 

1. f i SOS/SOSr, 

2. There exists a rational vector y G Q'" and the associated linear form Ljj G (]R[X]2d)* such 
that for any polynomials v,u E M.[X] with supp{v) G T and supp{u) G Qj^Tj ^e have 
mfv"^) <0 andLij{u^) > 0, 

where m is the number of elements in the set {X""^^ \ X°',X'^ G ^/,r}- 



4. Hilbert-Artin representation of positive semidefinite 
rational functions 

We generalize our method for solving the following problem: Given a rational function f / g G 
with g{X) ^ an integer e > and T C Terms [X; deg<(i], certify f /g ^ SOS/SOSr- 
Consider the following set 



[iyW,l^P]] 



g{X) ■ mTerms[X;deg<d]^W^'^^"iTerms[X;deg<d] = f {X) ■ niT^W^^^mr . 

lyW y 0, iy[2] y 0, Tr(l^[2]) ^ ^ ^ I J 



where d = e + ([deg(/) — deg{g)])/2. 

Proposition 4.1 We have f/g ^ SOS/SOS7- if and only if the set (12) is empty. 

Let 

Ti := {x"+^+^ I X^ e supp(^),X",X^ e Terms[X; deg<rf]} , 

T2 := {X"+/3+7 I x^ e supp(/), X", X^ G r} . 

We assume that Fi ^ r2, otherwise f/g ^ SOS/SOSr- The following analysis is similar to 
the one given in Section 2. 

The primal block SDP considered here has the same form as (7) but we use 

giX) ■ mTerms[X;deg<d]^W^f'WTerm.[X;dcg<d] = J^^G^"'^ * W^^'')^" (14) 



to define matrices G'"! 


Its dual problem is 




s* := inf s 

(j/,s)eM™+i 




s.t.M{y,s)yO, 


where 


M{y,s):= 


' Maigy) 

M,{{-f)y) + si 



(15) 



y '■= {ya)aeN" £ M.^'^d and m is the number of elements in the set Fi. M^^gy) and Me{{—f)y) 
are localizing moment matrices. Similar to Lemma 2.3, we have 

Lemma 4.2 There exists p = {y, s) such that Md{gy) >- and Me{—fy) + si >~ 0. 

Proof. Taking tja to be the one defined in the proof of Lemma 2.3, since g{X) is nonnegative, 
for any polynomial q{X) G M.[X] with supp(g) G Terms[X; deg<d], we have 

{vec{q)'^ , Md{gy)vec{q)) = / g{x)q{xffi{dx) 

g{x)q{xYh{x)dx 
> whenever g 7^ 0, 

10 



which imphes M^lgy) >- 0. Take 

then Me{-fy) + sJ ^ 0. D 

Based on the Farkas Lemma 2.4 and Lemma 4.2, similar to Theorem 2.5 and Theorem 
2.7, we have the following results. 

Theorem 4.3 Given a rational function f / g G Q(X) with g{X) >z and an integer e > 0, 
let d = e + ([deg(/) — deg(5f)])/2, then for any subset T C Terms[X; deg<e], the following 
are equivalent: 

1. f/g i SOS/SOSr, 

2. Fi ^ F2 in (13), or there exists a rational vector y = (y^) G Q'" such that Md{gy) ^ 
and Meify) -< in (15), 

where m is the number of elements in the set Fi. 

With a view towards linear forms in M[X], Theorem 4.3 is equivalent to 

Theorem 4.4 Given a rational function f/g G Q(X) with g{X) >z and an integer e > 0, 
let 'R[X]2d+dcg{g) '■= {p G ]R[X] I supp{p) G Terms[X; deg<2(i + deg(5f)]} where d = e + 
([deg(/)— deg(5')])/2, then for any subsetT C Terms[X; deg<e], the following are equivalent: 

1. f/g i SOS/SOSr, 

2. Fi ^ F2 in (13), or there exists a rational vector y G Q"* and the associated linear form 
L^ G (]R[X]2d+dcg(g))*; such that for any polynomials f,M G M.[X] with suppiy) G T and 
supp{u) G Terms[X; deg<d], we have Ljj^fv"^) < and L^i^gu"^) > 0, 

where m is the number of elements in the set Fi. 

5. Examples and Experiments 

Example 5.1 We prove that the well known Motzkin polynomial 

fix,, X2) = XfX2 + X^X^^ + 1 - 3X^X2 

is not an SOS. Set n = 2, e = and d = 3. By exploiting the sparsity, we have ^/,deg<o = 
{l,XiX2,XfX2,XiX|}. According to Corollary 3.3, m = 10 and we need to find a rational 
sequence y G Q^° or its associated linear form Ly G (MfXjg)* such that for any polynomial 
u G ]R[X] with supp(m) G ^/,dcg<05 we have Ly{v?) > and Ly{f) < 0. The certificate we 
obtained is 

y =(yo,o = yi,i = yi,2 = o, ^2,2 = 300, ^3^2 = ^2,3 = ^4,2 = ^3,3 = m^ = o)- 

Its associated linear form Ly satisfies 

L^(u2) = 300n^ 1 > 0, 
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and 

Ly{f) = -3 X 300 = -900 < 0, 

which imphes / can not be written as an SOS. 

Example 5.2 In [Kaltofen, Yang, and Zhi 2009], the monotone column permanent (MCP) 
conjecture has been proven for dimension 4 via certifying polynomials Pi,i, ^1,2, Pi,3, P2,2, P2,3, P3,3 
of degree 8 in 8 variables to be positive semidefinite, see [Kaltofen, Yang, and Zhi 2009] for 
the explicit forms of these polynomials. Among them, the polynomials Pi, 1,^3,3 are perfect 
squares. Applying the hybrid symbolic-numeric algorithm in [Kaltofen, Li, Yang, and Zhi 
2012], they proved that the polynomial pi,3 can be written as an SOS and the polynomials 
Pi,2,P2,2,P2,3 can be written as an SOS divided by weighted sums of squares of variables. 
We certify that all polynomials Pi,2,P2,2,P2,3 can not be written as an SOS via finding the 
corresponding certificates y G Q™ and the associated linear forms L^ G (]R[X]8)*. By ex- 
ploiting the sparsity, for pi^2, the matrices W, M{y,s) in (7), (8) are of dimension 24 x 24 
and m = 189. For p2,2, W and M{y, s) are of dimension 29 x 29 and m = 255. For p2,3, W 
and M{y, s) are of dimension 39 x 39 and m = 372. 

Example 5.3 This example comes from the even symmetric sextics in [Choi, Lam, and 
Reznick 1987]. Let 

n 

for integers k, < k < n — 1, we define polynomials fn,k by 

/„,o = -nMn,6 + in + l)Mn,2MnA - M^^, 

and 

fn,k = (A;' + k)Mn,6 - {2k + 1)M„,2M„,4 + M^_2, 1 < A; < n - 1. 

Some interesting results about these polynomials have been given in [Choi, Lam, and Reznick 

1987]. 

Proposition 5.1 For n > 3, 

(1) all fn,k, < k < n — 1, are positive semidefinite polynomials; 

(2) the polynomials fn,o and fn,i are SOS; 

(3) the polynomials fn,2, ■ ■ ■ , fn,n-i are not SOS; 

(4) M3,2 ■ /3,2 is an SOS [Robinson 1973]; Ml^ ■ /4,2 is an SOS; 

(5) for n > 4, Mn,2 ■ fn,n~i is an SOS. 

For n > 4 and 2 < i < n — 2, we wish to know whether Mn,2 ■ fn,i is an SOS. We have the 
following results. 



Ex.5. 3.1 For n = 4, we can certify that the polynomial /4^2 ^ S0S/S0Sdeg<2- By exploiting 

the sparsity, in (7), W^^^ has dimension 55 x 55 and W^'^'^ has dimension 5x5. We 
have m = 369 in (8). 
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Ex.5. 3. 2 For n = 5, we can certify the following: 



/5,2 ^ S0S/S0Sdeg<2 and /s.s ^ SOS/SOS 



dcg<4- 



By exploiting the sparsity, for /5^2, l/^t^l, 14^'^] have dimension 105 x 105,6 x 6, 
respectively and m = 1036. For /s^s, "H^I^l, "H^P] j^g^^g dimension 231 x 231, 21 x 21, 
respectively and m = 2751. 

Ex. 5. 3. 3 For n = 6, we can certify the following: 



/6,2 ^ S0S/S0Sdcg<2 and fe,3, /6,4 ^ S0S/S0Sdeg<4- 



By exploiting the sparsity, for /6,2, P^t^l,!^'^' have dimension 182 x 182,7 x 7, 
respectively and m = 2541. For /g^a and /6,4, Vr'^l,!^'^^ have dimension 434 x 
434, 28 X 28, respectively and m = 7546. 

Example 5.4 Consider the polynomial f{Xi,X2) = Xf + X| - 2X1X2 = (Xi - Xs)^ Its 
minimum is 0. However, for any small perturbation e > 0, the polynomial /e(Xi,X2) = 
(1 - e^)Xf + X| - 2X1X2 is not an SOS. Indeed, fe{C, C) = -e^C^ which implies that the 
infimum of /^ is —00. Hence / is an z/Z-posec? polynomial [Button, Kaltofen, and Zhi 2010]. 
For e = 10~^, . . . , 10~^, we can use Matlab SDP solver SeDuMi in Step III in Algorithm 3.1 
to certify that fe is not SOS. But for e < 10^^, Step III does not work out and we are not able 
to obtain a rational solution at which Sk < 0. If we use the command findsos in SOSTOOLS 
[Prajna, Papachristodoulou, and Parrilo 2002], it outputs a wrong SOS decomposition. Our 
method implemented in SDPTools in Maple can give exact certificates for /^ being not an 
SOS for e = 10^^ or smaller! Take e = 10^^ for instance. By exploiting the sparsity, we have 
^/e,dcg<o = {^15^2}- Setting Digits = 45 in Maple, the certificate we obtained is 

46635362642387337096986 . 53470001073377890290267 



y = Z/2,0 = -,^.o1^o^101ooor^^ror1^^r '^1.1 



1731626131338905851065 " ' 1985404333861113854675 
19926414238854847715525\ 



1/0,2 



739891310902398542446 / ' 
For any u G M.[X] with supp(M) G ^/^,dcg<0! we have 

2n _ 46635362642387337096986 ^ 19926414238854847715525 ^ 

^^^ ' ~ 1731626131338905851065 ^^'° ^ 739891310902398542446 "°'^ 
53470001073377890290267 

+ 2 X Ui nUn 1 

1985404333861113854675 ' ' 
, //46635362642387337096986 19926414238854847715525\ ^ 

> \2Ui nUn 1 X 

- ' '" "' ' \^V 1731626131338905851065 739891310902398542446 J 
53470001073377890290267\ 



1985404333861113854675 



) 



>0. 



However, 



9999999999999999 46635362642387337096986 19926414238854847715525 

^^■'^^ "10000000000000000 ^ 1731626131338905851065 ^ 739891310902398542446 
53470001073377890290267 

^ 1985404333861113854675 ^ ' 
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which imphes /^ is not SOS. 

Example 5.5 In this example, we consider some rational functions. 

Ex. 5. 5.1 For Motzkin polynomial in Example 5.1, we can certify that 

XtXl + XlXl + l-?,XlXl ^ gog/gog^^^^^_ 
-^l + i 

Ex. 5. 5. 2 For the even symmetric sextics in Example 5.3, we can certify that 

^,...,%^^S0S/S0Sdeg<2, n = 4,5,6. 

Those constitute our largest certificates. 
The correctness of the above result is guaranteed by the following proposition. 

Proposition 5.2 Let f/g G ]R(X) he a multivariate rational function where f,g & ]R[X] 
with GCB{f,g) = 1. ///, -/ ^ SOS then f/g i S0S/S0Sdeg<deg(9) ■ 



Proof. Assume the contrary, namely that 



9 EU^Axy 



u,{X),v^{X) e R[X], deg(t;,) < deg{g)/2. (16) 



Thus the right-side of (16) constitutes the reduced fraction f/g, which means g{X) = 
c^- Vj{Xy for some non-zero constant c G M, making f/c = Ei=iUi{XY, a contradic- 
tion, n 
Furthermore, for the polynomials in Example 5.3, we can compute the certificates for the 
following result: 

^"•' ^ S0S/S0Sdeg<4, n = 4,5,6 and Al ^ S0S/S0Sdeg<6. (17) 



We have no generalization of Proposition 5.2 to ±/ ^ SOS/SOSdcg<2e5 and the impossibilities 
(17) may hint of new unknown properties of the even symmetric sextics in [Choi, Lam, and 
Reznick 1987]. 
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